Center domains and their phenomenological consequences 
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We argue that the domain structure of deconfined QCD matter, which can be inferred from the 
properties of the Polyakov loop, can simultaneously explain the two most prominent experimentally 
verified features of the quark-gluon plasma, namely its large opacity as well as its near ideal fluid 
properties. 
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One of the major achievements of the experimental 
program at the Relativistic Heavy Ion Collider (RHIC) 
is the creation of a novel state of hot and dense QCD mat- 
ter dubbed the strongly interacting Quark-Gluon Plasma 
(sQGP) [IH5]. The characterization of the properties 
of the sQGP is based on three major discoveries: (1) 
the measurement of strong elliptic flow and the success 
of relativistic viscous hydrodynamics with a very small 
value of the shear viscosity r\ to entropy density s ratio 
close to the conjectured quantum lower bound [7] that 
argues for the near perfect liquid nature of the sQGP 
and a very short thermalization time of less than 1 fm/c, 
(2) the measurement of very strong suppression of high- 
momentum particles and rapid redistribution of the jet 
energy into the whole solid angle, which is indicative of 
the large opacity of the produced matter, and (3) the ob- 
served constituent quark number scaling law for the ellip- 
tic flow of identified hadrons as predicted by the parton 
recombination model [5HTT] that provides the most direct 
evidence for the formation of hadrons from a deconfined 
system of interacting patrons. All three discoveries have 
by now been confirmed in measurements of Pb+Pb col- 
lisions at the CERN Large Hadron Collider (LHC) [T2]. 

We note that observations (1) and (2) do not provide 
direct evidence for quark deconfincmcnt in the sQGP 
state, but argue in favor of a strongly coupled nature of 
the state. On the other hand, (3) provides the most di- 
rect evidence to date for quark deconfinement, but sheds 
little light on the dynamical properties of the sQGP. At 
present, there exists no rigorous theoretical approach in 
the framework of quantum chromodynamics (QCD) that 
could explain all of the above features. Most calculations 
of the jet energy-loss are based on perturbative QCD and 
are only sensitive to the gluon content of the matter, but 
do not distinguish between different microscopic struc- 
tures related to quark confinement or chiral symmetry 
breaking. Lattice calculations are not yet able to give 
reliable results for transport properties of the sQGP that 
can provide clues about its structure. Holographic ap- 
proaches to strongly coupled supersymmetric gauge the- 
ories have been used to model dynamical properties of 
the sQGP, but it remains unclear how well these gauge 
plasmas mirror the physics of the QGP. 

So far, very little attention has been paid to the non- 
perturbative structure of the gauge field configurations 
in the quark-gluon plasma and how it may affect the 



features observed in the experiments at RHIC and LHC. 
In particular, our focus here is on the Polyakov loop. 
We shall argue that the sQGP produced in relativistic 
heavy ion collisions has a domain structure based on the 
different minima of the Polyakov loop potential in the 
deconfined phase, and that this domain structure can be 
instrumental for generating the large opacity and small 
value of rj/s seen in the sQGP. 

Two characteristic features distinguish the quark- 
gluon plasma from the hadron phase: restored chiral 
symmetry in the light quark sector and color decon- 
finement. The order parameter which characterizes the 
confinement-deconfincment transition in the pure SU(3) 
gauge theory is the Polyakov loop L. It is defined as 



L(x) = — tr P exp 
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where P is an ordering operator with regard to the imag- 
inary time r, A 4 = iA° = iA 0a \ a /2 with the Gell-Mann 
matrices A a (a = 1, • • • , 8), T the temperature, and x 
the three-dimensional spatial position. In the pure gauge 
theory, the Polyakov loop and the free energy of a heavy 
static quark, Fq(T), are related by 



F Q (x,T) = -T\n\(L(x, T))|, 



(2) 



where (•) stands for the thermal average. (L(x, T)) van- 
ishes in the confined phase. In the deconfined phase it 
takes the value of one of the three elements of Z3, the 
center of SU(3): exp(i2^?r/3) with v = 0, 1, 2. Thus, the 
confined phase is Z3 symmetric, while the Zj, symme- 
try is spontaneously broken in the hot deconfined phase. 
Lattice calculations in the deconfined phase are usually 
carried out around one of these three states. Here we ex- 
plore how the center symmetry affects the properties of 
the hot matter created in relativistic heavy ion collisions. 

The fundamental Polyakov loop |l]), = L, gov- 
erns the interaction of a static quark with the thermal 
gauge field. The interaction of particles carrying color 
charge in the adjoint representation, such as gluons, is 
described by the adjoint Polyakov loop. There is good 
evidence from lattice simulations [121 [T5] that the expec- 
tation value of the adjoint Polyakov loop is given 
by a power of (L^) equal to the ratio of the Casimir 
operators of the adjoint and fundamental representation 
of SU(3): (L^) w (L( 3 )) 9 / 4 . This means that gluons 
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interact even more strongly with the center domain walls 
than quarks. These results are for the pure gauge theory. 
In the early stage of heavy ion collisions, as we discuss 
in the following, the system is gluon dominated and the 
effect of dynamical quarks is expected to be suppressed. 

It is now widely believed that the quark-gluon plasma 
is formed through a non-equilibrated precursor state 
called glasma [15] . The glasma is a longitudinally slowly 
varying, classical gauge field configuration whose struc- 
ture in the transverse direction is characterized by the 
saturation scale Q s , with Q s sa 1.5 GeV at RHIC 
and Q s m 2 GeV at LHC [16]. Through thermaliza- 
tion this gauge field materializes into the quark-gluon 
plasma. The values of Polyakov loop in the plasma will 
locally cluster around the one of the preferred values 
exp(z2^7r/3) where the potential has a minimum. Since 
the original transverse correlation length of the gauge 
fields is of order Qj 1 , causality dictates that it can be 
at most of order rth at the moment of thcrmalization. 
In other words, the distribution of Polyakov loop values 
must assume a domain structure on the transverse plane 
with a typical transverse domain size Rd constrained by 
Qj 1 < Rd < T t h. Each domain is separated from other 
domains by potential walls, where the Polyakov loop val- 
ues interpolate between the different Z 3 values. We call 
these domains center domains |17l I18j . 
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FIG. 1: In the deconfined phase of QCD, center domains are 
characterized by one of the three values exp(i2z/7r/3) of the 
expectation value of the Polyakov loop (L) corresponding to 
the minima of the Polyakov loop potential. In the boundary 
wall between two domains, (L) changes from one minimum to 
another. The figure shows two types of schematic trajectories. 

There are two possible types of trajectories for the 
Polyakov loop between two Z 3 minima. Either the expec- 
tation value of the Polyakov loop remains of unit mod- 
ulus, but its phase changes gradually by 2-7r/3 inside the 
domain wall (solid path in Fig. [IJ, or the modulus of the 
Polyakov loop becomes smaller than unity and abruptly 
changes its phase near the origin (dashed path in Fig[T]). 
In order to determine which case is realized, we need the 
information of the Polyakov loop potential in the decon- 
fined phase. Its form has been empirically determined by 
a combination of analytical considerations and by fitting 
the lattice results for the pressure as a function of temper- 



ature. A commonly used form in the so-called Polyakov 
loop-Nambu-Jona-Lasinio (PNJL) model is [19] 

U((L)) = -bT[54e- a / T \(L)\ 2 +lnP((L),(tf))] . (3) 

with P(z,z) = 1 — 6|z| 2 - 3|z| 4 + 4(z 3 + z 3 ) arising 
from the SU(3) Haar measure and numerical constants 
a = 0.664 GeV and b = 0.0075 GeV 3 . T denotes the 
temperature. Contour plots of the Polyakov loop poten- 
tial ([3]) at three different temperatures below, near and 
above the deconfinement temperature (100, 300, and 500 
MeV) are shown in Fig. [2] Since the potential U((L)) 
is deduced from lattice simulations, its empirical form is 
most reliable in the vicinity of the potential minima, i.e. 
(L) = exp(i2^7r/3), but it is not well known far away 
from these minima. The Polyakov loop potential is thus 
not sufficiently accurately determined to reliably answer 
the question with certainty, which trajectory (L) takes 
inside the domain wall. 
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FIG. 2: Contour plots of the SU(3) Polyakov loop potential 
§ at T = 100 MeV (left), T = 300 MeV (middle), and 
T = 500 MeV (right). The transition from the confined to the 
deconfined phase is clearly visible by a shift of the minimum 
from (L) = to (L) = exp(z2i/7r/3). 

Assuming the widely adopted potential form (J3|, the 
Polyakov loop prefers to take values around zero inside 
the domain walls between deconfined regions, because 
trajectories corresponding to the dashed path in Fig. [T] 
involve the lowest potential barrier. Since the expecta- 
tion value of the Polyakov loop in the confined phase 
vanishes, center domains are thus separated from neigh- 
boring ones by walls that are characterized by gauge field 
configurations similar to those in the confined phase of 
the gauge theory. This insight forms the basis of our 
following discussion. 

We now proceed to consider the possible phenomeno- 
logical consequences of the presence of the center domains 
in the quark-gluon plasma produced in heavy ion colli- 
sions. First, the walls act as a potential barrier for par- 
tons with momentum smaller than the confining scale. 
This can be seen by invoking the following argument. 
The free energy of a single static quark is related to the 
expectation value of the Polyakov loop by the well known 
expression A vanishing average Polyakov loop thus 
corresponds to an infinite energy of an isolated heavy 
quark. For light quarks, this argument does not imply 
confinement because of the possibility of pair creation, 
but it implies that it is energetically unfavorable for any 
quark to propagate through a region with (L) w 0. 
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Thus, partons with thermal momenta cannot classi- 
cally penetrate the walls but instead are reflected on 
them. Consequently, the mean free path Af of partons 
becomes of order i?d/2. Inserting this result into the ki- 
netic theory formula for the shear viscosity, 



V 



rnpXf 



(4) 



where n denotes the particle density and p is the thermal 
momentum, one obtains 



TRd, 



(5) 



for p 3T, and s/n « 4 valid for massless particles. 
Because of the small domain size, the shear viscosity is 
highly suppressed. This mechanism of lowering the shear 
viscosity resembles that of the anomalous viscosity [2UI 
12 lj , but the mechanism responsible for a small mean free 
path is different. 

For quarks with higher momenta, which can cross the 
boundary between domains, the system of domain walls 
acts like the combination of a frequency collimator and 
an irregular undulator. Because they serve as a reflective 
barrier to soft gluons, the domain walls constitute a very 
effective mechanism for the "frequency collimation" effect 
|22j . which has been conjectured to be responsible for the 
large di-jet asymmetries observed in Pb+Pb collisions at 
the LHC [23], [24] . The stripped soft gluon component 
of the nascent jet is quickly restored, because the walls 
alternately decelerate and then re-accelerate fast moving 
partons. This can be seen by expanding the expression of 
the average of the Polyakov loop and making a Gaussian 
approximation: 
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Inside the domain walls, where the Polyakov loop ap- 
proximately vanishes, the gauge field A fluctuates with 
large amplitude, and partons thus feel sudden changes 
of the zeroth component of the external vector potential 
through the coupling in the covariant derivative. This 
leads to strong and uncorrelated radiation of gluons ev- 
ery time a parton crosses one of the walls, which are dis- 
tributed randomly with the density of the order of R^ 3 - 
The center domains thus provide a novel, nonpertur- 
bative mechanism for jet quenching, which is different 
from other mechanisms in two aspects. First, it distin- 
guishes between the QCD phases: jet quenching by cen- 
ter domains is caused by the confining phase composing 
the domain walls. In the conventional mechanisms, jet 
quenching is not affected by which phase the matter is in, 
but only by how large the gluon density is. Second, the 
gluon radiation is very different at low and high frequen- 
cies: At frequencies higher than the height of the domain 



wall barrier, u> > lo c , gluons propagate as quasi-particles, 
but radiation of gluons is enhanced by partially coherent 
emission from the crossing of many domain walls. At low 
frequencies, radiated gluons are strongly scattered by the 
domain walls and propagate diffusively, not ballistically. 

To obtain an estimate of the energy loss in a single 
wall crossing we simply assume that all gluons below a 
critical frequency uj c are scattered away from the parent 
parton. The energy loss is then given by 
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(7) 



where dN g /duj denotes the frequency spectrum of glu- 
ons accompanying the hard parton. If just the Lienart- 
Wiechert field of the hard parton is restored after a wall 
crossing, the gluon spectrum is given in the Weizsacker- 
Williams approximation by: 



dN^ 
dio 



C 2 a s 



ln(w/wo) 9{uj — loq) 



(8) 



where we have introduced an infrared cut-off loq ~ 0(T) 
accounting for the finite size of the domains, which lim- 
its the ability of the parton to radiate very soft gluons. 
Ci denotes the Casimir for the color charge of the hard 
parton. The energy loss per unit length is then 



dE 
dx 



C 2 a s 
irR d 



{w c ln(w c /wo) - (w c - w )} 



(9) 



where Rd is the average domain size. Estimating the 
cut-off frequency as w c ss 1 — 2 GeV, ujo ~ 0.4 GeV, 
and Rd ~ 0.5 fm, we find values of dE / dx in the range 
(0.2 - 1) C 2 a s GeV/fm. 

In addition to this energy loss associated with energy 
collimation, for gluon frequencies above w c , radiation is 
enhanced by the interaction of the energetic parton with 
several domain walls. Making use of the undulator anal- 
ogy, the characteristic frequency of radiation emitted by 
a massive parton traversing the system of domain walls 
is given by [15] 



UJ L = AlTJ 2 /R d , 



(10) 



where 7 = E/m is the Lorentz factor of the parton. Obvi- 
ously, these nonperturbative effects can lead to enormous 
energy loss of moderately heavy quarks, such as charm 
quarks, at energies a few times their rest mass. 

Another consequence of the center domain scenario is 
that the momenta of the emitted soft gluons are imme- 
diately randomized by the reflection on the walls. This 
property can explain the almost complete redistribution 
of the quenched energy over the whole solid angle ob- 
served in Pb+Pb collisions at LHC [26]. Figure [3] depicts 
two partons propagating through a small volume of QGP 
matter: a hard parton radiates (soft) gluons while cross- 
ing domain walls whereas soft partons may reflect off the 
domain walls. 

The time evolution of the center domains has been 
studied in a different context [27) . Gupta et al. assumed 
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FIG. 3: Schematic representation of the domain structure of 
a small volume of QGP matter with two patrons traversing 
it. A fast parton crossing the walls of a center domain will 
radiate (soft) gluons that may reflect from the domain walls, 
leading to large energy loss and rapid isotropization of the jet 
energy. Soft partons may reflect multiple times off domain 
walls while propagating through the QGP. 

that center domains are created by the tunneling of the 
Polyakov loop expectation values in the early stage of 
the quark-gluon plasma formation. The expected size 
of the domains is then much larger than Qj 1 , and it is 
assumed to possess the boost-invariant structure for sim- 
plicity. According to their numerical simulation, domains 
merge if their initial values of (L) lie in the same mini- 
mum, while domains do not merge if their initial values 
of (L) correspond to different minima. 

Applying the results of ref. |27| to our scenario, we con- 
clude that both the size and number of the initially cre- 



ated domains does not change much during the time evo- 
lution of the quark-gluon plasma. The mean free path of 
partons does not increase substantially beyond the range 
of Qj 1 estimated above, and the correlation length the 
gauge fields is limited by this scale, while the system is 
strongly interacting. This implies that hydrodynamics 
remains applicable throughout the expansion until the 
deconfinement transition. 

When the temperature approaches the pseudo-critical 
temperature T c , the potential difference between (L) = 
and (L) = exp(i2^7r/3) becomes smaller, the walls be- 
come wider, and the deconfined domains become smaller. 
Eventually, deconfined domains fragment into individual 
hadrons, and the walls become the nonperturbative QCD 
vacuum. The center domains thus describe the evolution 
of the quark-gluon plasma consistently from its birth to 
hadronization and help to explain the strongly coupled 
nature of the quark-gluon plasma including all the ob- 
served properties of the quark-gluon plasma from its hy- 
drodynamical behavior to jet quenching. 

In summary, we have argued that the center domains 
are an important facet of the evolution of the quark-gluon 
plasma from its birth up to hadronization. They natu- 
rally explain the strongly coupled nature of the quark- 
gluon plasma including its major observed properties 
from its nearly ideal hydrodynamical behavior to strong 
jet quenching. 

B. M. and S. A. B. acknowledge support by U.S. De- 
partment of Energy grants DE-FG02-05ER41367 and de- 
sc0005396. M. A. was supported by Grants- in- Aid for 
Scientific Research 23540307 from JSPS of Japan. 



I. Arsene et al. (BRAHMS), Nucl. Phys. A 757, f (2005). [15 

K. Adcox et al. (PHENIX), Nucl. Phys. A757, 184 

(2005). [16 

B. B. Back et al., Nucl. Phys. A757, 28 (2005). 

J. Adams et al. (STAR), Nucl. Phys. A757, 102 (2005). [17 

M. Gyulassy and L. McLerran, Nucl. Phys. A750, 30 [18 

(2005). 

B. Miiller and J. L. Nagle, Ann. Rev. Nucl. Part. Sci. 56, [19 
93 (2006). [20 
P. Kovtun, D. T. Son, and A. O. Starinets, Phys. Rev. 
Lett. 94, 111601 (2005). [21 
R. J. Fries, B. Miiller, C. Nonaka, and S. A. Bass, Phys. 
Rev. C68, 044902 (2003). [22 
V. Greco, C. M. Ko, and P. Levai, Phys. Rev. C68, 
034904 (2003). [23 
D. Molnar and S. A. Voloshin, Phys. Rev. Lett. 91, 
092301 (2003). [24 
R. J. Fries, V. Greco, and P. Sorensen, 
Ann.Rev.Nucl.Part.Sci. 58, 177 (2008). [25 
B. Miiller, J. Schukraft, and B. Wyslouch (2012), 
1202.3233. [26 
S. Gupta, K. Huebner, and O. Kaczmarek, Nucl. Phys. 
A785, 278 (2007). [27 
A. Mykkanen, M. Panero, and K. Rummukainen, JHEP 
1205, 069 (2012), 1202.2762. 



T. Lappi and L. McLerran, Nucl. Phys. A772, 200 
(2006). 

J. L. Albacete, A. Dumitru, and Y. Nara, 

J.Phys.Conf.Ser. 316, 012011 (2011). 

R. D. Pisarski, Phys.Rev. D62, 111501 (2000). 

A. Dumitru and R. D. Pisarski, Phys.Lett. B504, 282 

(2001). 

K. Fukushima, Phys.Lett. B591, 277 (2004). 

M. Asakawa, S. A. Bass, and B. Miiller, Phys. Rev. Lett. 

96, 252301 (2006). 

M. Asakawa, S. A. Bass, and B. Miiller, 
Prog.Theor.Phys. 116, 725 (2007). 

J. Casalderrey-Solana, J. G. Milhano, and U. A. Wiede- 
mann, J.Phys.G G38, 035006 (2011). 
G. Aad et al. (Atlas Collaboration), Phys.Rev. Lett. 105, 
252303 (2010). 

S. Chatrchyan et al. (CMS Collaboration), Phys.Rev. 
C84, 024906 (2011). 

C. Brau, Free electron lasers, Academic Press, San Diego 
(1990). 

M. B. Tonjes (CMS Collaboration), J.Phys.G G38, 
124084 (2011). 

U. S. Gupta, R. K. Mohapatra, A. M. Srivastava, and 
V. K. Tiwari, Phys.Rev. D82, 074020 (2010). 



